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Abstract 

The influence of Coulomb interaction on the electron transport through molecular 
wires is studied in the regime of incoherent tunneling. In the limit of strong Coulomb 
repulsion, the current for spinless electrons is determined. It is shown that the 
voltage profile along the wire crucially influences the dependence of the current on 
the wire length. Upon inclusion of the spin degree of freedom one finds a blocking 
eff'ect which depends both on the interaction strength and on the population of the 
wire. For finite Coulomb interaction, the temperature dependence of the blocking 
is studied and it is shown that several regimes with different blocking strength may 
exist. 
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1 Introduction 



In the last few years, considerable progress towards an electronics based on 
single molecules has been made [1]. Nanotubes [2,3], fullerenes [4,5], and or- 
ganic molecules [6-9] have been contacted and current-voltage characteristics 
were measured. These experimental successes are accompanied by a vast body 
of theoretical studies [10] which mostly concentrate on coherent transport 
through the molecules. However, relaxation processes on the molecule may 
lead to incoherent transport [11-13] and some of the experimental data may 
be interpreted in this sense [5,8]. 
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Fig. 1. The molecular wire is modeled by identical wire units represented by circles 
which are joined by identical molecular bridges "M" . Two contacts "C" provide the 
connection to the left and right electrodes denoted by "L" and "R" , respectively. 

Recently, Petrov et al. [12,13] have applied a nonlinear quantum kinetic sta- 
tistical many-body approach in mean-field approximation to incoherent trans- 
port through short molecular wires in the limit of infinitely strong Coulomb 
interaction. Stimulated by this work, we discuss the charge transport through 
a single wire in the incoherent regime. Assuming that relaxation processes 
on the wire units are much faster than the time scale for tunneling between 
different units, a description in terms of rate equations for the electronic state 
of the wire is then adequate. 

The main focus of this paper will be on the influence of Coulomb interaction 
on the electronic transport. After introducing in Section 2 a generic model for 
the molecular wire and discussing the voltage profile across the wire, the limit 
of strong Coulomb repulsion will be considered in Section 3. This allows us to 
restrict the calculation to the lowest unoccupied molecular orbitals (LUMO) 
and to neglect excitations to higher levels. In particular, this imphes that 
temperature is sufficiently low so that the Coulomb interaction exceeds the 
thermal energy. In Section 4 we will include the spin degree of freedom for the 
electrons which had been neglected in Section 3. In this context the blocking of 
spin channels due to Coulomb interaction will be discussed in detail. Finally, in 
Section 5 we relax the requirement of strong Coulomb interaction and study 
the different regimes arising from the interplay between Coulomb repulsion 
and thermal fluctuations. 



2 Model for a molecular wire 

In this paper, we will focus on generic properties of incoherent charge trans- 
port through a molecular wire. The model for the wire will consequently be 
kept as simple as possible: a linear array of N identical wire units numbered 
consecutively from left to right is joined by — 1 identical molecular bridges 
denoted by "M" . In order to apply a voltage and to pass a current through 
the wire, the setup shown in Fig. 1 is used where two contacts "C" connect 
the wire to the electrodes marked by "L" and "R" . 

Except for Section 5, we will conflne the discussion to the case of strong 
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Fig. 2. The energies of the LUMOs of the wire units are depicted as dashed Hnes 
at energy £ above the electrochemical potential in the left electrode. An external 
voltage V shifts the levels to the positions marked by the full lines. The energy 
differences A^Jl, ^Em, and Ai?R are defined in Eqs. (2)-(4). ryL, ??M, and ?7r denote 
the fractions of the total applied voltage which drop across the left contact, the 
wire, and the right contact, respectively. 

Coulomb interaction. This prevents the wire from being charged with more 
than one electron. Neglecting spin for the moment, the wire will then always 
be found in one of iV + 1 states: either the wire is in its ground state and all 
wire units are neutral, or exactly one of the units is charged by an additional 
electron. Wc will denote these states by "0" for a neutral wire or the number 
of the unit i = 1, . . . , N which carries an extra charge. 

Since we are interested in the current across the molecule beyond linear re- 
sponse, we also need to specify the voltage drop along the wire. In princi- 
ple, this requires the solution of the coupled Schrodinger and Poisson equa- 
tions of the electrode- molecule-electrode system [14], which in itself presents 
a formidable task that we will not try to tackle. Instead, we will use the sim- 
ple model proposed in Refs. [7,15], which views the two electrodes as parallel 
plates forming a capacitor with the molecule in between acting like a dielectric. 

Within this picture one obtains a potential profile where the applied voltage 
V drops partially by rjiy and rj-^V at the left and right contact, respectively, 
and decreases linearly by rji^V along the wire. The condition 77l + ^71^ + = 1 
ensures that we have accounted for the full voltage drop. This situation is 
shown in Fig. 2 where we have fixed the electrochemical potential in the left 
electrode at zero energy. The dashed lines at energy £ indicate the position of 
the LUMOs in the absence of an applied voltage. The actual levels represented 
by the full hnes at energies 



are shifted according to the local voltage, thereby displaying the voltage pro- 
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Fig. 3. The positions of the energy levels axe shown for models A (left) and B (right) 
of the voltage drop across the system. 

The current through the wire is governed by the energy differences between 
neighboring sites 



and between the electrochemical potentials in the electrodes and the adjacent 
wire units 



Particularly at very low temperatures the current depends strongly on the 
energy differences (2)-(4). If one of these energy differences is negative, the 
absence of thermally activated processes will lead to a suppression of the 
current. This situation is depicted in Fig. 2 where AEi^ < 0. If, on the other 
hand, the energy levels of the molecular sites lie between the electrochemical 
potentials in the two electrodes, a current may flow even at zero temperature. 
Apart from Section 5, we will always have the latter situation in mind when 
discussing zero temperature results. 

In the following, the electrochemical potential —eV in the right electrode is 
assumed to lie lower in energy, thereby favoring transport from left to right, 
which we will denote as forward direction. Furthermore, we specialize the 
voltage proflle to the two limiting cases depicted in Fig. 3. Within model A 
the voltage is constant along the molecule and drops by equal amounts at the 
two contacts (?7l = ''7r = 1/2, t^m = 0). In contrast, in model B the entire 
voltage drops linearly along the molecule (?7l = '7r = 0, t/m = !)• 



3 Current for spinless electrons in the presence of strong Coulomb 
interaction 



/S.Em = Ei — Ei+i = Tju 
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Fig. 4. Transitions between the N + 1 states of the wire occur either by tunneling 
between adjacent wire units or between the end units of the wire and one of the 

two electrodes. 

3.1 Master equation and tunneling rates 



The incoherent electron transport through an individual wire shall be de- 
scribed as sequential tunneling events of an electron between the electrodes 
and the end units of the wire or between adjacent wire units. The dynamics 
is then modeled in terms of a Markovian master equation. Fig. 4 shows the 
possible transitions between the + 1 different wire states. The arrows over 
the corresponding tunneling rates F indicate the direction of transport for a 
geometry depicted in Fig. 1, while the subscripts L, M, and R correspond to 
the left contact, a molecular bridge, and the right contact, respectively. 

In order to determine the average current through the wire, we will need the 
probabilities Pi{t) to find the system in one of the states i = 0, . . . ,N. Since 
these constitute the only possible states, the probabilities should satisfy the 
normalization condition 

= (5) 

1=0 

The dynamical evolution of the probabilities due to the incoherent transitions 
depicted in Fig. 4 is determined by the master equation 

Po{t) = -(f L + f R)po(t) + f LPl(t) + f Rpjv(t) (6) 

Pi{t) = -(f M + f L)pi(t) + r^poit) + rMP2(t) (7) 

Pi{t) ^-{fM + fu)Pi{t) + fuPi-i{t) + fuPi+i{t) {i^2,...,N -1) (8) 

PN(t) = -(Fr + rM)Piv(i) + f RPo(i) + ruPN-i(t) . (9) 
From a solution of these equations, the average currents through the left and 
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hit) = e(rLPo - TlPi), /rW = e{TnPN - TrPo), (10) 
and through the molecular bridge between units i and i + 1 

lM,i{'t) = e(rMPi - FMPi+l) (11) 

can be obtained. We will be particularly interested in the stationary case where 
the probabilities pj become time-independent. Then, all currents (10) and (11) 
will be equal, i.e. 

I = h = I^ = Iu,i , (12) 
thereby ensuring current conservation. 

As can be seen from Eqs. (6)-(9), the stationary probabilities pi depend only 
on ratios of transition rates. We therefore introduce the "backward-forward" 
ratios 

Tm Tl Tr 

7m = — , 7l = — , 7r = — , (13) 
Tm Tl Tr 

and the "branching" ratios 

c.=h, c.=h. (14) 

The latter two relate molecular properties to contact properties and therefore 
depend on the microscopic details. The ratios (13), on the other hand, combine 
backward and forward rates at the same contact or bridge. In order to ensure 
thermodynamic equilibrium they are given by Boltzmann factors at inverse 
temperature (3 = l/k-QT 

7M = e-/^^^M ^ ^ e-^^^^ , 7R = e"^^^- , (15) 

where the energy differences have been introduced in Eqs. (2)-(4). Since all 
energy differences have to sum up to the difference in electrochemical potential 
between left and right electrode, i.e. 

AEi^ + {N- 1)AEm + AE^ = eV , (16) 

the "backward-forward" ratios obey the relation 

7l7^-Sr = e-^^^ . (17) 
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3.2 Stationary solution 



In the following, we will restrict ourselves to the stationary case. Together 
with the ratios (13) and (14), the master equation (6)-(9) then turns into 

-(CR + 7RCL)P0 + 7LCRPl + CLPiV = (18) 

-(Cl + 7l)Pi +Po + 7mClP2 = (19) 
-(1 + 7M)Pi + 7MPi+i = (i = 2, . . . , AT - 1) (20) 

-(1 + ^mCr)pn + 7rPo + CrPn-1 = . (21) 

Out of these + 1 equations only N are linearly independent, so that the 
normalization condition (5) is needed to determine the N +1 stationary prob- 
abilities Pi- 

It is convenient to start with the "molecular" part (20), which takes on the 
form of a two-step recursion relation and can thus readily be solved in terms 
of the two probabilities pi and pn leading to 

1 _ ^N-i _ N-i 

Pi = -. N^Pi TT^Pn ^ = l,...,iV . 22 

1 - 7m 1 - 7m 

Inserting this result into Eqs. (19) and (21) and making use of the normaliza- 
tion (5), one finds for the probability for a neutral wire 



Po = 



JV 1 "I — 1 

(7^ - 1)(Cl7r + Cr) + N{e-(^^y -1) + ^(1 - 7l7r) 
^ 7l(7^-^ - 1) + (7M - 1)(Cl + Cr7l7^-') 



(23) 



while the probability for an electron on site i may be expressed in terms of po 
as 



_ e-^^^ -1 + [(7M - 1)(Cl7r + Cr) + 1 - 7l7r] 7m"' 

Pi ~ I N-l i\ , / I A N-l\ Po ■ V^^J 

7l(7m -1) + (7m-1)(Cl + Cr7l7m ) 



3.3 Stationary current 



The stationary solution (23) and (24) allows to calculate from any of the 
expressions (10) and (11) the average current (12). The resulting expression 
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can be cast into the form 



1 _ Q-PeV 
Tl Tm 

where we have introduced the abbreviations 

CL = 1 + ^ ^4l^, (26) 



7l 1 - 1/7: 



M 



cm = 7lV-^ + (7l7r-1) ' , (27) 

1 - 7m (1 - 7m)^ 1 - 7m 

CR=^ + (28) 

7r 1 - 7m 

The current changes its sign when the apphed voltage is reversed and at the 
same time forward and backward rates as well as left and right contacts 
are interchanged. This can be verified by noting that under these changes 
the coefficients (26)-(28) are transformed according to cl/Fl g^^^cr/Fr, 
cr/Fr ^ e/^^^ cl/F L, and cm/Fm ^ e^^^^ cm/Fm- 

Further insight into the result (25)-(28) can be gained by discussing several 
special cases. An important limit, which will be particularly relevant in the 
discussion of the following sections, is a wire consisting only of one unit and 
no bridges. Then, the current simplifies to 



FlFr-FlFr ^29) 



Fl + Fl + Fr+Fr 



containing forward and backward contributions. For positive energy differences 
A£^L and AE^, the backward rates vanish at zero temperature and one is left 
with the familiar result 



e (30) 
Fl + Fr 



for an sequential incoherent process with rates Fl and Fr. 

For more than one wire unit it is instructive to compare the two models for the 
voltage profile introduced in Section 2. We start by considering model A where 
the voltage drops only across the contacts. As a consequence, the forward and 
backward rates on the molecule are equal and we have 7m = 1- In this case. 
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the coefficients (26)-(28) become 



(31) 



Cl = 1 + 



CM = 7L(iV-l) + (l + e-^^^) 



N{N-1) 
2 



(32) 



Q-/3eV 

CR = + N 

7r 



(33) 



which is consistent with Eqs. (29) and (30) for = 1. For long wires, the 
current scales like l/N"^ with the number of wire units. Since forward and 
backward tunneling rates on the wire are equal, the electron has equal proba- 
bility for tunneling in forward and backward direction. The length dependence 
of the current can therefore be understood in terms of a random walk where 
the average time to diffuse across a wire of length N is proportional to N'^. 

For model B, where the voltage drops entirely along the wire, the same ar- 
gument holds provided the thermal energy is much larger than the energy 
difference between adjacent sites, ksT » AE'm- At low temperatures, how- 
ever, the situation differs qualitatively since the backward rates are negligibly 
small. In order to obtain a finite current even in the limit of zero temperature, 
we assume that AE'l, A£^r, > 0. At zero temperature, Eqs. 25)-(28) then yield 



which reduces to (30) for N — 1. For a long wire, the current now decreases 
linearly with the number of units or, equivalently, with the length of the wire. 
In this sense, for model B the wire resembles an ohmic resistor. 

The comparison of models A and B demonstrates that the dependence of 
the current on the wire length at low temperatures may provide qualitative 
information about the voltage drop across the molecule. In any case, however, 
incoherent transport leads to an algebraic length dependence of the current 
in clear contrast to the exponential dependence characteristic for coherent 
transport [10]. 



4 Blocking of spin channels by strong Coulomb interaction 

So far, we have ignored the spin degree of freedom of the electrons. In this 
section, we will take the spin into account but, due to the strong Coulomb 
interaction, there will be not more than one electron present on the wire. 
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Fig. 5. The situation depicted in Fig. 4 is generalized to include the electron spin. 
The two rings corresponding to the two spin directions are joined via the neutral 
state "0" of the molecule. 

Furthermore, we will assume that the two spin directions are equivalent and 
that no spin flips occur. 

Inclusion of the electron spin increases the number of possible wire states to 
2N + 1. The wire can now either be neutral (with probability Poit)) or there 
can be an extra electron with spin up or down at site i — 1, . . . ,N (with 
probability Pi,i{t) and Pi,i{t), respectively). The possible transitions between 
the different states arc depicted in Fig. 5. We emphasize that, like in Section 3, 
all rates are to be understood as rates per spin direction. In the presence of a 
magnetic field, the two spin directions would no longer be equivalent and the 
rates would have to be distinguished by a spin index. 

A comparison of Figs. 4 and 5 shows that the two molecular parts are indepen- 
dent of the spin direction and correspond to the molecular part in the spinless 
case. The only difference lies in the coupling to the neutral state "0". This 
motivates to introduce the total probability pi{t) = Pi,i{t) +Pi,i{t) for an elec- 
tron on site i irrespective of its spin. The resulting master equation including 
the normalization condition is of the same form as Eqs. (5) and (6)-(9) in the 
spinless case. The only difference consists in a factor of 2 by which the rates 
Fl and describing the charging of the wire have to be multiplied. For these 
rates, the spin degree of freedom leads to an additional process allowing to 
occupy the wire. This is not the case for an electron leaving the wire. In this 
case, the spin is fixed and the corresponding rates Fl and Fr are not doubled. 
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The appearance of a factor of 2 may also be understood in terms of the density 
of states in the electrodes which enters the tunneling rates. For tunneling of 
electrons onto the wire the number of initial states is increased by the spin 
degree of freedom thereby leading to an effective doubhng of the density of 
states. On the other hand, for tunneling into the electrodes the density of final 
states is not increased since the spin is determined by the electron leaving the 
wire. 

Taking into account the factor of 2 in the tunneling rates that populate the 
wire, one finds for the average current 

/ = e(2rLPo-fLPi) . (35) 

Here, the probability pi implicitly contains a factor of 2 since the electron con- 
tributing to the backward current may carry either spin up or down. Making 

the required modifications of the rates Fl and Fr, the current (25)-(28) turns 
into 



1 - e-f^'V 
Ch I Qvi ^ cr 

Fl Fm Fr, 



(36) 



where the coefficients now are given by 



1 , e-^^y 1 - l/7f 

2 7l 1 - I/tm 



CL = J + ^ ^^T^, (37) 



2 1 - 7m (1 - ImV 1 - 7m 

Cr = -7, ^ ■ (39) 

27r 1 - 7m 

We remark, that due to the mutual blocking of the electrons one does not 
obtain an overall factor of 2 in the current relative to the spinless case (25)- 
(28). Instead, the current shows a more intricate infiuence of the spin degree 
of freedom. 



The mechanism of blocking is particularly transparent if the wire consists of 
only one site, i.e. N = 1. We restrict ourselves to the case of zero temperature 
where, as in the derivation of Eq. (30), all backward rates vanish. Then the 
current (36) simplifies to read 

/ = 2e J^^l . (40) 
2Fl + Fr 

This result has to be compared with the case where the two spin directions lead 
to independent transport channels so that the expression (30) is multiplied by 
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Fig. 6. The current suppression (42) due to Coulomb interaction is plotted as a 
function of the ratio rj of the tunneling rates in the left and right contacts. 



a factor of 2 yielding the current 



In = 2e: 



Li R 



Tl + Tr 



(41) 



The effect of blocking due to the Coulomb interaction is quantified by the 
ratio of the current (40) in the presence of Coulomb repulsion and the current 
(41) in the absence of interaction, i.e. 



l_ _ 1 + 7] 

Io~ l + 2ri' 
The blocking is completely determined by the ratio 



(42) 



V 



Tr 



(43) 



of the rates through the left and right contacts and varies between 1 for Fl <^ 
Fr and 1/2 for Fl » Fr as shown in Fig. 6. 

This dependence on 77 can be understood in terms of the average population of 
the wire which here is given by pi. From the stationary solution of the master 
equation for the case with spin and N = 1 one obtains 



Pi 



2r] 



(44) 



l + 2r) 

The blocking effect on the current (42) may therefore be directly related to 
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the average population of the wire by 



For small rj, or equivalently Fl ^ Fr, the wire is depleted almost immediately 
after it has been populated. The average population of the wire therefore is 
very small and the passage of another electron is almost never blocked. In the 
opposite case of large r], or Fl ^ Fr, the right contact represents a bottleneck 
which leads to a sizable average population of the wire and therefore to the 
blocking of transport. 

This result clearly shows, that inclusion of the spin degree of freedom does not 
necessarily lead to a doubling of the current as would be the case in the absence 
of Coulomb interaction on the wire. Apart from the interaction strength, the 
average population of the wire plays a decisive role and blocking is strongest 
when the right contact increases the population by acting as a bottleneck. 



5 Temperature dependence of blocking 



The calculation of the current-voltage characteristics presented so far was 
based on the assumption that at any given time there is at most one extra 
electron on the wire. Although it is straightforward to generalize the rate 
equations (6)-(9) to more than one electron by extending the state space, in 
general one has to solve the rate equations numerically. 

Some insight into the validity of the limit of strong Coulomb interaction can 
be gained by allowing a second electron on the wire. Then, an analytical 
treatment is still possible if we restrict the length of the wire to one site, i. e. 
N — 1. We will assume that this very site may be occupied by at most two 
electrons of opposite spin. Then, there exist four different wire states, |0), |t), 
I J,), and III), which correspond to no electron, one electron with spin up or 
down, and two electrons, respectively. The energies of these four states are 
schematically shown in Fig. 7 where we allow for an interaction energy U if 
two electrons are present on the wire. By changing both, the temperature as 
well as the one-electron energy E relative to the left electrochemical potential, 
we may choose the levels which are relevant for the transport through the 
wire. 

The discussion in the last section has shown that the effect of blocking is 
most pronounced when the tunneling rates through the right contact are much 
smaller than through the left one. We therefore shall concentrate on this limit. 
Then, the master equation for the probabilities p„ of being in a state with n — 
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Fig. 7. Many-particle states for one site with n = 0, 1, and 2 electrons. A is the 
one-particle energy of the states |t) and \i) and U is the Coulomb interaction 
energy between the electrons occupying the site. 

0, 1, or 2 electrons depends only on the rates through the left contact. These 
tunneling rates may be identified by the change of the number of electrons. 
For example, r2^i refers to the backward rate through the left contact with 
initially two electrons on the molecule. With this notation the master equation 
reads 

Po{t) = -2ro^ipo(t) + ri^oPi(t) (46) 
Pi{t) = -(Fi^o + ^i^2)pi{t) + 2ro^iPo(0 + 2r2_iP2(i) (47) 
P2(t) = -2r2^ip2(i) + ri_2Pi(t) (48) 

with the normalization condition 

Po{t) + Pi{t) + P2{t) = I . (49) 

The rates Fq^i and F2_>i are each multiplied by a factor of 2, thereby account- 
ing for the two spin directions possible in these processes. The probability pi 
of finding one electron of arbitrary spin on the wire implicitly contains a factor 
of 2 in the same sense as was discussed after Eq. (35). 

If the probabihties are calculated from the master equation (46)-(48), thereby 
neglecting tunneling through the right contact, the current has to be evaluated 
at that contact. We will assume that the electrochemical potential in the right 
electrode is so low that the forward rates Fr through the right contact are 
approximately independent of the occupation of the molecular site and that 
backward rates may be neglected. Then, the stationary current is obtained as 

/ = erR(pi + 2p2) , (50) 

where the factor of 2 accounts for the two spin directions available when one 
of two electrons is leaving the molecule. 

Solving the master equation (46)-(48) in the stationary limit, we obtain from 
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(50) for the current 



ri^2 
i + f^ 

/ = eVn p . (51) 

2ro_>i 2r2_*i 

In order to study the effect of blocking, we have to compare this current with 
the current Iq in the absence of Coulomb interaction. In this case ri_»2 — ^o-^i 
and r2_»i = ri_>o and we find 

lo = eFR . (52) 

2ro^i 2ri^o 

This result still depends on temperature and on the energy £ of the one- 
electron level relative to the left electrochemical potential. 

The ratio I/Iq depends only on rate ratios which, as in Section 3.1, may be 
expressed in terms of Boltzmann factors. For single occupation of the molecule, 
we have 

^ = e-^' ■ (53) 
Double occupancy costs the additional interaction energy U, so that 

^ = e-^(^+^) . (54) 

r2->i 

The blocking due to Coulomb interaction is then described by 

L- 1 - e-^^ 



Figure 8 depicts this result as a function of kBT/\£\ for U/\£\ = 2 and 2000. 
The full lines correspond to the case, where the one-electron energy lies below 
the left electrochemical potential, i.e. £ < 0. The dashed line refers to the 
opposite case with £ > 0. As a function of temperature, we can distinguish 
three different regimes. 

In the limit of zero temperature, one finds for £ < 0, and U > —£ an inter- 
action induced reduction of the current by a factor 1/2. This is the blocking 
discussed already in Section 4. For £ + U < 0, double occupancy is possible 
even in the zero temperature limit and no blocking occurs. As can be seen 
from the dashed fines in Fig. 8, blocking is also absent for £^ > 0. In this case. 
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Fig. 8. The current (51) through the wire normahzed with respect to the current (52) 
in the absence of Coulomb interaction is plotted as a function of the temperature 
T in units of \£\ for U/\£\ = 2 and 2000. For each value of U/\£\, a full line and a 
dashed line are shown, corresponding to £ < and £ > 0, respectively. 

the current as well as the occupation of the molecule approach zero exponen- 
tially for low temperatures. In analogy to the reasoning presented in Section 4, 
blocking thus becomes ineffective. 

For U ^ \S\, a second regime appears at intermediate temperatures \S\ <C 
ksT <^ U. For U/\£\ — 2000, Fig. 8 indeed shows a clear plateau at a blocking 
factor of 2/3. This can be attributed to the fact that for fcgT ^ \S\ the 
backward rate Fi^o is of the order of the forward rate Fq^i. This reduces 
the population of the molecule and therefore blocking is less effective than at 
lower temperatures. 

The reduction by a factor 2/3 can be understood even without recourse to 
Eq. (55). For equal forward and backward rates, the probabilities for an empty 
wire or a wire occupied with one electron of given spin are equal, leading to 
probabilities po = 1/3, pi = 2/3 and P2 = 0- From Eq. (50) one therefore finds 

/ = (2/3)eFR,. In the absence of Coulomb interaction, double occupancy is 
allowed and one finds all four possible states with equal probability leading to 
Po = 1/4, pi = 1/2 and p2 = 1/4. This yields the current Iq — bFr, and thus 
the blocking factor reads I/Iq — 2/3. 

Finally, in the high-temperature regime, where ksT 3> the thermal 

energy is much bigger than the interaction energy U so that the latter becomes 
irrelevant. Therefore, the Coulomb interaction can no longer lead to blocking 
of the electronic transport. 
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6 Conclusions 



The influence of Coulomb interaction on incoherent electronic transport 
through a molecular wire has been found to be multifaceted. In the regime of 
strong Coulomb interaction, the current of spinless electrons shows an alge- 
braic dependence on the wire length characteristic for incoherent transport. 
It was found, however, that this length dependence varies as a function of the 
voltage proflle across the system. At low temperatures a, behavior indi- 
cates the absence of a voltage drop along the molecule while in the presence 
of a voltage drop a. 1/N behavior is expected. 

Inclusion of the spin degree of freedom in the absence of Coulomb interaction 
leads to an increase of the current by a factor of 2. Coulomb interaction, on 

the other hand, should lead to a reduction of the current. It was found that 
this blocking effect not only depends on the interaction strength but also on 
the average population of the wire and thus on the tunneling rates (cf. Fig. 6). 

For finite Coulomb interaction, the amount of blocking also depends on tem- 
perature. Between the limits of suppression of the current by a factor of 2 
at zero temperature and the absence of blocking at high temperatures an 
intermediate regime may exist. There, Coulomb interaction prohibits double 
occupancy. Nevertheless the current is only reduced by a factor of 2/3 (cf. 
Fig. 8). This is a consequence of the competition between backward and for- 
ward rates, which equally populates the two one-electron levels. 

As a final remark, we mention that even though we have restricted ourselves 
to electronic transport, the approach presented here is also applicable to hole 

transport. Depending on the situation, holes may be the dominant charge 
carriers; the results presented above may easily be adapted to this case. On 
the other hand, if both, holes and electrons, contribute to the transport, an 
extension similar to the one presented in Section 4, where spin is included, 
becomes necessary. 



Acknowledgements 

The authors would like to acknowledge useful discussions with I. Goychuk. 
This work has been supported by the Deutsche Forschungsgemeinschaft 
through Sonderforschungsbereich 486, project AlO. 

References 



17 



[I] C. Joachim, J.K. Gimzewski, A. Aviram, Nature 408 (2000) 541. 

[2] T.W. Ebbesen, H.J. Lezec, H. Hiura, J.W. Bennett, H.F. Ghaemi, T. Thio, 
Nature 382 (1996) 54. 

[3] S.J. Tans, M.H. Devoret, H.J. Dai, A. Thess, R.E. Smalley, L.J. Geerligs, 
C. Dekker, Nature 386 (1997) 474. 

[4] C. Joachim, J.K. Gimzewski, R.R. Schhttler, C. Chavy, Phys. Rev. Lett. 74 
(1995) 2102. 

[5] D. Porath, Y. Levi, M. Tarabiah, O. Millo, Phys. Rev. B 56 (1997) 9829. 

[6] M.A. Reed, C. Zhou, C.J. Muller, T.P. Burgin, J.M. Tour, Science 278 (1997) 
252. 

[7] S. Datta, W. Tian, S. Hong, R. Reifenberger, J.I. Henderson, C.P. Kubiak, 
Phys. Rev. Lett. 79 (1997) 2530. 

[8] C. Kergueris, J.-P. Bourgoin, S. Palacin, D. Esteve, C. Urbina, M. Magoga, 
C. Joachim, Phys. Rev. B 59 (1999) 12505. 

[9] J. Reichert, R. Ochs, D. Beckmann, H.B. Weber, M. Mayor, H. v. Lohneysen, 
|cond-mat/0106219 . 

[10] for a review see e.g. A. Nitzan, Ann. Rev. Phys. Chem. 52 (2001) 681. 

[II] D. Segal, A. Nitzan, W.B. Davis, M.R. Wasielewski, M.A. Ratner, J. Phys. 
Chem. B 104 (2000) 3817. 

[12] E.G. Petrov, P. Hanggi, Phys. Rev. Lett. 86 (2001) 2862. 

[13] E.G. Petrov, V. May, P. Hanggi, this issue. 

[14] V. Mujica, A.E. Roitberg, M. Ratner, J. Chem. Phys. 112 (2000) 6834. 

[15] W. Tian, S. Datta, S. Hong, R. Reifenberger, J.I. Henderson, C.P. Kubiak, J. 
Chem. Phys. 109 (1998) 2874. 



18 



